Numerical investigation of photon creation in a three-dimensional resonantly 

vibrating cavity: TE-modes 
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The creation of TE-mode photons in a three-dimensional perfectly conducting cavity with one 
resonantly vibrating wall is studied numerically. We show that the creation of TE-mode photons 
in a rectangular cavity is related to the production of massive scalar particles on a time-dependent 
interval. The equations of motion are solved numerically which allows to take into account the 
intermode coupling. We compare the numerical results with analytical predictions and discuss the 
effects of the intermode coupling in detail. The numerical simulations reveal that photon creation 
in a three-dimensional resonantly vibrating cavity can be maximized by arranging the size of the 
cavity such that certain conditions are realized. In particular, the creation of TE-mode photons 
in the lowest frequency mode (1,1,1) is most efficient in a non-cubic cavity where the size of the 
non-dynamical dimensions is roughly 11 times larger than the size of the dynamical dimension. We 
discuss this effect and its relation to the intermode coupling in detail. 

PACS numbers: 03.65.-w, 03.70.-fk, 12.20.Ds, 42.50.Lc 
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I. INTRODUCTION 

In 1948 Casimir Q predicted an attractive force be- 
tween two perfectly conducting plates (ideal mirrors). 
This so-called Casimir effect |S H, ^ IS caused by the 
change of the zero point energy of the quantized elec- 
tromagnetic field in the presence of boundaries has been 
verified experimentally with high accuracy |E HIE IS 01- 
The existence of the Casimir force ^3 acting on macro- 
scopic boundaries confirms the reality of quantum vac- 
uum fluctuations and their potential influence even on 
macroscopic scales. 

Besides the change of the zero point energy of the quan- 
tum vacuum provoked by static boundary conditions a 
second and even more fascinating feature of the quantum 
vacuum appears when considering dynamical, i.e time- 
dependent boundary conditions. The quantum vacuum 
responds to time- varying boundaries with the creation of 
real particles (photons) out of virtual quantum vacuum 
fluctuations. This effect, usually referred to as dynamical 
or non-stationary Casimir effect '12^ , has gained growing 
interest during recent years. 

A scenario of particular interest are so-called vibrat- 
ing cavities where the distance between two parallel 
mirrors changes periodically in time. The possibility of 
resonance effects between the mechanical motion of the 
mirror and the quantum vacuum leading to an even ex- 
ponential growth of the particle occupation numbers for 
the resonance modes makes this configuration the most 
promising candidate for an experimental verification of 
the dynamical Casimir effect. 

For a one-dimensional vibratin g ca vit y this effect has 

been studied in numerous work Ji|JiHjilJil El El 

mmEllillilllllillll Elllllsaiar showing 



'Electronic address: Marcus. Ruser@physics.unige.ch 



that the total energy inside a resonantly vibrating cavity 
increases exponentially in time. 

The more realistic case of a three-dimensional cavity is 
studied in H H S S El H El • The important 
difference between one- and higher-dimensional cavities 
is that the frequency spectrum in only one spatial dimen- 
sion is equidistant while it is in general non-equidistant 
for more spatial dimensions. An equidistant spectrum 
yields strong intermode coupling whereas in case of a non- 
equidistant spectrum only a few or even no modes may 
be coupled allowing for exponential photon creation in a 
resonantly vibrating three-dimensional cavity [lHl33 . l37j . 
Without intermode coupling the equations of motion for 
the field modes reduce to harmonic oscillators with time- 
dependent frequency. Particle creation can then be inves- 
tigated by using an approach based on Schrodinger scat- 
tering theory [40j. Even though for higher-dimensional 
cavities the problem can be reduced to a single harmonic 
oscillator in some special cases Ell the intermode cou- 
pling cannot be neglected in g eneral [ssL Is^ . (Se e also 
the discussion of the work MlJj in section IX of [l5l|.) 

Field quantization inside cavities with non-perfect 
boundary conditions has been studied in, e.g., 0, 
and corrections due to finite temperature effects are 
treated in m0,El|- The interaction between the quan- 
tum vacuum and the (classical) dynamics of the cavity 
has been investigated in [2l |H IH El and an approach 
to the dynamical Casimir effect based on stationary walls 
but time-dependent conductivity properties is discussed 
in _4||. 

The electromagnetic field inside a dynamical cavity can 
be decomposed into components corresponding to the 
electric field parallel or perpendicular to the moving mir- 
ror. It is then possible to introduce vector potentials for 
each polarization, transverse electric (TE) and transverse 
magnetic (TM) ,,3^ ,5flj, ^Slj . The equations of motion for 
TE-modes in a dynamical rectangular cavity are equiv- 
alent to the equations of motion for a scalar field with 
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(time-dependent) Dirichlet boundary conditions |33l l37| . 
More complicated boundary conditions, so-called gen- 
eralized Neumann boundary conditions, emerge when 
studying TM-modes 's^', . In most of the works cited 
above only TE-polarizations arc treated. For recent work 
dealing also with TM-polarizations see [s^ IH^ ] . 

The aim of the present work is to study photon cre- 
ation in a vibrating three-dimensional cavity fully numer- 
ically taking the intermode coupling into account. We 
show that the equations of motion describing the dynam- 
ics of the transverse electric modes (TE) in a dynamical 
rectangular cavity correspond to the equations of motion 
for a massive scalar field on a time-dependent interval 
(one-dimensional cavity). Thereby the wave number of 
the TE-modes associated with the non-dynamical cav- 
ity dimensions is identified with the mass of the scalar 
field. Creation of TE-polarized photons can then be stud- 
ied with a formalism presented and tested for a mass- 
less scalar field in a one-dimensional cavity in 53-, 54]. 
Even though the method of [sM Is^ is valid for a vari- 
ety of boundary conditions it is not directly applicable to 
generalized Neumann boundary conditions which involve 
a time-derivative appearing when studying TM-modes. 
For other recent numerical work see also [ssl Is^ Is^ . 

The paper is organized as follows. In section II we 
present the equations of motion for TE-modes in a three- 
dimensional rectangular cavity and show that they corre- 
spond to the equations of motion for a massive scalar field 
in a one-dimensional cavity. The formalism for study- 
ing the dynamical Casimir effect for a massive scalar 
field on a time-dependent interval numerically is reviewed 
in section III. Some analytical results obtained for TE- 
mode photons are summarized in section IV. We present 
and interpret the numerical results in section V and 
discuss their consequences for photon creation in three- 
dimensional vibrating cavities in section VI. We conclude 
in section VII and discuss some details about the numer- 
ics in the appendix. 



II. EQUATIONS OF MOTION FOR TE-MODES 
IN A RECTANGULAR DYNAMICAL CAVITY 

The dynamics of the transverse electric modes (TE- 
modes) inside a rectangular ideal (i.e. perfectly con- 
ducting) cavity of dimensions {(0, Ix), (0, ly){0, h)} is de- 
scribed by the wave (Klein-Gordon) equation ^ 



[a2-A]<i>(t,x) = o 



(1) 



with the massless scalar field $(t,x) subject to Dirichlet 
boundary conditions at all walls of the cavity [3^ . l3^ . 

The x-diniension of the cavity is assumed to be dynam- 
ical with the right wall following a prescribed trajectory 



l(t) = lx{t). At any moment in time the field can be 
expanded as 



$(<,x) =^g„(i)0„(<,x) 



with canonical variables qn{t) and functions 



(2) 




(3) 



ensuring Dirichlet boundary conditions at the positions 
of the cavity walls |33j . The functions 4>n{ii x) form an or- 
thonormal and complete set of instantaneous eigenfunc- 
tions of the Laplacian A with time-dependent eigenvalues 



(4) 




Each field mode is labeled by three integers nx,ny,nz = 
1, 2, ... for which we use the abbreviation n = (ux, ny,nz). 

Inserting the expansion Q into the field equation 
multiplying it with (prnit, x) and integrating over the spa- 
tial dimensions leads to the equation of motion for the 
canonical variables qn{t) 33]: 

qn{t) + nl{t)qn{t) + 2Y,M^n{t)qm{t) (5) 

m 

+ [^«>n(i) - Nnmit)] qm{t) = 0. 

m 

The time-dependent coupling matrices Afnm(^) and 
Nnin{t) are given by [13 

lit) 



(6) 

m r (-1)"-+"^ ^^j^^^j^^^^ if nx ^ fux 

l{t) \ ii Ux ^ fUx 



and 



k 



(7) 



^ We are using units with h - 



During the dynamics of the mirror the time evolution of 
a field mode n may be coupled to (even infinite many) 
other modes m via the time-dependent coupling ma- 
trix Afnm(i)- In Eq. ((SJ for a given mode {nx,ny,nz) 
the coupling matrix (O yields couplings of g(n^,n„.n.) to 
g(m,,n„,«,) and q(m^,ny.n,), i-s. Only summations over irix 
appear. Modes with different quantum numbers in the 
y— and z— directions are not coupled and the quantum 
numbers corresponding to the non-dynamical dimensions 
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enter the equations of motion only globally. Therefore we 
can identify q^it) = (?(„^,„^,„^) (i) and 



with n = Tlx and the wave number 



m 



+ fcy (8) 




(9) 



associated with the non-dynamical cavity dimensions. 

Because all summations over m ~ {mx,'my,mz) in- 
volving the coupling matrix (jnj reduce to summations 
over a single quantum number to, Eq. Q is equiva- 
lent to the differential equation describing a real massive 
scalar field on a time-dependent interval [0,Z(<)] (one- 
dimensional cavity) when fcy is identified with the mass 
of the field [s^l- As pointed out in ^37^ the number of 
created TE-mode photons equals the number of created 
Dirichlet scalar particles in a three-dimensional cavity. 
Consequently, the number of TE-mode photons created 
in a three-dimensional cavity equals the number of scalar 
particles of "mass" k\\ created in a one-dimensional cav- 
ity [0,^(t)]. Photon production in TE-modes in a three- 
dimensional cavity can therefore be studied numerically 
with the formalism presented in |^ . 



III. THE FORMALISM 

Quantization is achieved by replacing the set of classi- 
cal canonical variables {qn,Pm\ with the corresponding 
operators {qmPm} and demanding the usual equal-time 
commutation relations. Furthermore, the Heisenberg pic- 
ture is adopted from now on. The relation between the 
canonical variable 5„ and the canonical momentum is 
given by p„ ^„ + Y,m ImMmn- Assuming that the cav- 
ity is at rest for times t < the coupling matrix vanishes 
and equation jSJ reduces to the equation of a harmonic 
oscillator with constant frequency il^ = 51„(i < 0). Con- 
sequently 



quit < 0) 



with frequency 



1 



a„e 



-in?, t 



(10) 



(11) 



where Ig = /(O) and we have introduced the dimensionless 
"mass parameter" M — lok\\. The time-independent an- 
nihilation and creation operators an,aJi associated with 
the particle notion for t < are subject to the commu- 
tation relations 



0, \an,al 



(12) 



The initial vacuum state \0,t < 0) is defined by 

a„|0,i<0) =OVn. (13) 

When the cavity dynamics is switched on at f = and 
the wall follows the prescribed trajectory l{t) field modes 
are coupled due to the non-vanishing coupling matrix 
Mnm- To account for the coupling, the operator q„ may 
be expanded as [s^ 



1 



a„e(™)(i)+a]„e(™)*(i)l (14) 



with complex functions eir'\t) satisfying Eq. If the 
motion ceases and the wall is at rest again for t > ti the 
operator g„(t > ti) takes the form ^ 



qn{t > ti) 



(15) 

with = fln{t > ti) and annihilation and creation 
operators A„, AJ^ corresponding to the particle notion for 
t > ti. The final vacuum state |0, t > ti) is defined by 



An\0,t>ti) =OVn. 



(16) 



The initial state particle operators a„, ajj are linked to the 
final state particle operators An , A\ by the Bogoliubov 
transformation 



(17) 



and the number of particles (photons) created in a mode 
n during the motion of the wall is given by the expecta- 
tion value of the number operator A'lAn associated with 
the particle notion for t > ti with respect to the initial 
vacuum state |0,t < 0): 



iV«(ii) = (0,t<0|Ati„|0,t<0) 



(18) 



The total number of created particles as the sum of 
Nn(ti) over all quantum numbers n 

Nih) = J2 ^"(^l) = E E \^^nnih)\' (19) 



is in general ill-defined and requires appropriate regular- 
ization. This can be done most easily by introducing an 
explicit frequency cut-off which also simulates non-ideal 
boundary conditions for high frequency modes 14J. Asa 



^ Here /(ti) = li is assumed to be arbitrary. For an oscillating 
cavity, however, it is natural to consider times ti after which the 
dynamical wall has returned to its initial position. 
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matter of course, such a frequency cut-off has to be used 
in the numerical simulations. 

In order to calculate Bmniti) we introduce auxiliary 
fimctions and i}n^\t) via [H3 | 



with 



k 



(20) 



ryi'")(i) = 6M(i)-^ 

(21) 

Using the second order differential equation (jSJ for 
(i) it is easily shown that those functions satisfy the 
follow ing: system of coupled first-order differential equa- 
tions 



E 

k 



(22) 



^i"^ (i) = (^)d'"^ (i) - a+„ (t)77i") {t)_ 

- E[^n.WerW + c;,(t)r,i")(t)] (23) 



with 



a„„(i) 



1 ± 



r!„(i) 



9P 



M„fc(t)±^Mfe„(t) 



For the coupling matrix (jSJ one finds in particular 



lit) 



kr 



l{ty ' n^-k^ 



It 



(24) 
(25) 

(26) 



if n 7^ fc and c^„(i) = 0. The advantage of this sys- 
tem of first-order differential equations relies on the fact 
that, besides the time-dependent frequency ri„(t), only 
the coupling matrix Mnk enters but neither its square 
Nnk nor its time derivative M„fc. 

By matching Eq. (fTljl with Eq. (fT5|l for q„(t) and the 
corresponding expressions for Pn{t) aX t ~ ti one finds 
the relations jS^] 



A,t(0 - 2 



1 ± 



(29) 



Demanding that the field be in its vacuum state 
|0, i < 0) as long as the mirror is at rest implies Amn{^) — 
Smn and BmniO) = 0. Accordingly the initial conditions 
for ^i™^ (t) and 77!™^ (t) read 



^i™)(0) = 2<5_, ry(™)(0) = 0. 



(30) 



By means of Eq. H28() the number of created massive 
scalar particles, or equivalently the number of created 
TE-mode photons, at time t = t\ can now be calculated 
by solving the system of differential equations formed by 
Eqs. (|22|l and H23|l numerically using standard numerics. 
For this we truncate the infinite sums by introducing a 
cut-off quantum number fcmax to make the system of dif- 
ferential equations suitable for numerical treatment. The 
system is evolved up to a final time tmax and the particle 
number 1)18(1 is calculated for several times in between, 
i.e. we interprete ti as a continuous variable such that 
the particle number 1(181) becomes a continuous function 
of time Consequently, the stability of the numerical re- 
sults has to be guaranteed which means that for the low- 
est modes n the numerical values for Nn (t) remain prac- 
tically unchanged under variation of fcmax- More details 
regarding the numerics are collected in the appendix. 



IV. KNOWN ANALYTICAL RESULTS 

In what follows, we consider the periodic trajectory 

l{t) = lo[l + esm{iot)] , e< 1, (31) 

for which it was found in 33] that two modes I and k are 
coupled whenever one of the conditions given by 



(32) 



is satisfied ^. In a resonantly vibrating cavity uj = 217" 
with not one of those conditions fulfilled the number of 
TE-mode photons created in the resonant mode n in- 
creases exponentially in time jSSj: 



Nn{t) = sinh^(n7„ et) with 7„ = — ^ ( ^ 



(33) 



Bmn{tl) 




(27) 



(28) 



^ Potential problems inherent in this procedure like the appearance 
of discontinuities in the velocity of the mirror motion occurring 
when calculating the particle number for times t for which lit) 7^ 
are discussed in jMl in detail. We come back to this in section 
V. 

^ Here and in the following we have translated the results for three- 
dimensional cavities to the case of massive scalar particles ac- 
cording to Eq. ISl . 
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By means of multiple scale analysis the authors of |33 
also studied the resonance case lu — 2f2° with two cou- 
pled modes n and k satisfying 



31^ 



(34) 



For the particular case n — I and k — 5 analytical ex- 
pressions for the number of TE-mode photons are derived 
in [3^ . Given a mode n we can couple it to a particular 
mode k by tuning the mass M (or equivalently fcy) such 
that the condition H34(l is fulfilled. It is important to note 
that coupling between modes does occur even if Eq. H32|) 
is detuned, i.e. if Eq. H32() is satisfied by the frequencies 
fll and only approximately. The particular case of 
two modes n and k satisfying 



(3 + = nl 



(35) 



without additional couplings to higher modes was stud- 
ied in Q^. For sufficiently small k (i.e., k < e) the two 
modes n and k are still resonantly coupled and the num- 
ber of particles produced in both modes increasing expo- 
nentially with time. 



V. NUMERICAL RESULTS 

A. Preliminary remarks 

In [5^ IH^ I we have employed the same formalism to 
study the creation of massless scalar particles in a one- 
dimensional vibrating cavity numerically. In this case the 
numerical results agree with analytical predictions ob- 
tained under the assumption e <C 1 demonstrating the re- 
liability of the numerics. The extension to massive scalar 
fields is straightforward. We set Ig = 1, i.e. all physical 
quantities with dimensions are measured with respect to 
the length scale and dimensionless quantities are used 
throughout. The amplitude of the oscillations is fixed to 
e = 0.001 guaranteeing accordance between the numeri- 
cal results and the analytical predictions in the massless 
case mip. 

As mentioned before, we calculate the particle number 
for arbitrary times even though the analytical expres- 
sions we are comparing the numerical results with are 
valid only for times after which the dynamical wall has 
returned to its initial position. In 53] it is shown that (for 
a vibrating cavity) this leads to oscillations in the particle 
number which are of negligibly small amplitude when the 
amplitude of the cavity oscillations itself is small (e ^ 1). 
We will briefiy come back to this question later on. Fur- 
thermore, when calculating the particle number at times 
t for which l(t) 7^ the used particle definition requires 
a matching of the solutions to expressions corresponding 
to the static cavity l{t) = 0. This discontinuity in the ve- 
locity of the mirror trajectory may give rise to spurious 
contributions to the total particle number. The cut-off 
fcniax automatically ensures that the total particle num- 
ber remains finite because it automatically smoothes the 
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FIG. 1; (color online) Number of particles created in the res- 
onance mode n — 1 for mass parameters Af — 0.2, 0.7, 2, and 
3.5 in comparison with the analytical prediction 13311 . 



motion. For a more detailed discussion see |5J| where 
it is shown that the influence of such discontinuities is 
negligibly small in the case of a cavity vibrating with 
sufficiently small amplitudes (e — 0.001). (For a detailed 
discussion of how the initial discontinuity in the mirror 
motion H31|l affects the particle creation see also H^.) 
The numerical results which are presented and discussed 
in the following are practically not affected by the above 
mentioned effects. 



B. Main resonance uj = 

In Fig. n the number Ni{t) of particles created in 
the resonant mode A; = 1 is shown for masses M = 
0.2,0.7,2 and 3.5 and compared to the analytical pre- 
diction Eq. lO ^. For M = 0.7, 2 and 3.5 the numerical 
results are well described by Eq. (|33|) which is valid pro- 
vided that the resonant mode n = 1 is not coupled to 
other modes. In case of the mass M — 0.2 the numerical 
result for Ni{t) disagrees with the analytical prediction 
(|33|l . This will be discussed in the following in detail. 
Figure[21shows the corresponding particle spectra at time 
t = 6700. One infers that for M = 3.5, 2 and 0.7 the 
mode which becomes excited most is indeed the resonant 
mode n = 1. However, also higher modes become ex- 
cited but the corresponding particle numbers are several 
orders of magnitude smaller than the number of particles 
created in the resonant mode. For M — 0.7, for exam- 
ple, the mode A: = 3 is clearly excited. Figure O shows the 
number of particles created in the modes fc = 1, 2 and 3 
for the mass parameter M = 0.7 in detail. The difference 



^ In this section we use the general notion "particles" for massive 
scalar particles, or equivalently, TE-mode photons. Furthermore 
we call M the mass of the particle, having in mind that it corre- 
sponds to the wave number fen for TE-mode photons. 
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FIG. 2; (color online) Particle spectrum for different mass 
parameters M = 3.5,2,0.7 and M = 0.2 at time t = 6700 
corresponding to Fig. The spectra are shown for fcmax = 10 
(dots) and fcmax = 20 (squares) to demonstrate numerical 
stability. 
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FIG. 3: (color online) (a) Number of particles created in the 
modes fc = 1, 2 and 3 for the mass parameter M = 0.7 cor- 
responding to the spectra shown in Fig. |5| Part (b) shows 
N3{t) and part (c) N2{t) in each case for the two resolutions 
At = 0.01 (circles) and At = 0.005 (solid lines). The particle 
numbers calculated for times at which the mirror has returned 
to its initial position are accented by "+" and the background 
motion is shown for comparison as well. 



in the numerical values of iVi and N3 is so large that the 
contribution of A'3 to the total particle number is neg- 
ligible such that N ~ Ni. From Fig. |31 (a) one could 
conclude that N2, i.e. the number of particles created 
in the mode k = 2, behaves in the same way as N3 but 
shows superimposed oscillations. However, Figs. |31 (b) 
and (c) provide a more detailed view on the time evo- 
lution of Nk{t) for modes k = 2 and 3. In Fig. O (a) 
the resolution in which the numerical results are shown 
is not sufficient in order to resolve the details which are 
visible in Figs. 01 (b) and (c). These high resolution pic- 
tures reveal that N3 increases exponentially in time with 
oscillations superimposed on an average particle number 



whereas N2 itself oscillates strongly with an amplitude 
negligibly small compared to A^i . 

The small scale oscillations in the particle numbers are 
correlated with the periodic motion of the mirror which 
we have depicted in Figs. |21 (b) and (c) as well. But one 
infers that the particle numbers show oscillations also 
when the expectation value (|18|l is calculated only for 
times at which the mirror has returned to its initial posi- 
tion Iq. Therefore the oscillations cannot be traced back 
exclusively to the fact that Eq. H18|l is calculated for ar- 
bitrary times which may be considered as unphysical. 

The observation that also higher modes become excited 
(even though they are very much suppressed) is explained 
by the fact that two modes k and I are coupled even if 
Eq. (|32|) is not exactly satisfied by the two frequencies 
nl and fl'^. For M = 0.7 the equation 3fl^ = fl^ has no 
solution for integer k. Thus taking Eq. H32|l as an exact 
equation only the resonant mode should become excited 
and particle creation should take place in the mode n = 1 
exclusively. Inserting M ~ 0.7 one finds the solution 
k ~ 3.07 which is apparently close enough to the integer 
value fc = 3 to excite that mode. For smaller values of 
M the solution of 3ili — 51° approaches the value fc = 3 
and one has to expect that for sufficiently small values 
of M the mode coupling becomes again so strong that 
Eq. (|33|) does no longer describe the numerical results. 
This is the case for M = 0.2 yielding k ^ 3.005 for which 
a strong coupling between the modes n = 1 and k = 3 
occurs. Furthermore, from Eq. H32|) and the coupling of 
n = 1 and fc = 3 follows 2175 + = which has 
I ~ 5.004 as solution, i.e. the mode A: = 3 is coupled to 
the mode 1 = 5. In the same way the mode 5 is coupled to 
the mode 7. Thus interpreting Eq. (|22Jl as ~ jfi^ ± fi^j 
explains the numerically computed particle spectrum (cf. 
Fig. 121) which shows similar features as the spectrum 
obtained for the massless case (cf Fig. 4 (b) of ^53] ). One 
observes that also even modes become excited (like also 
for M = 0.7) which is not the case for M = [ssj. 
These modes are dragged by the strongly excited modes 
(odd modes) and rapidly oscillate (like N2 for M = 0.7) 
with an amplitude several orders of magnitude smaller 
compared to Ni, N3 and N^. In Fig. ^ we show the 
number of particles created in the modes k = 1 to 5 for 
M = 0.2 to illustrate the just-stated. As for M = 0.7 the 
number of particles created in the odd modes increases 
exponentially showing oscillations superimposed on an 
average particle number while the number of particles 
created in the even modes k = 2 and k = 4 consists of 
oscillations only with amplitudes much smaller compared 
to the number of particles created in the odd modes. As 
in Fig. O the strongly oscillating behavior of the particle 
numbers for even modes is visible in high time resolution 
only [part (c) of Fig.^. 

The fact that mode coupling occurs even if Eq. H32I) is 
not satisfied exactly is well known. We can rewrite the 
expression 311„ ~ fig to get (3 -I- K)fl^ = [Eq. 
As mentioned at the end of the former section it was 
shown for this case in Is^l that for sufficiently small k 
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FIG. 4: (color online) (a) Number of particles created in the 
modes fc = 1, 2, 3, 4 and 5 for the mass parameter AI = 0.2 cor- 
responding to the spectrum shown in Fig. |5| Part (b) shows 
N5{t) and part (c) N2(t) in each case for the two resolutions 
At = 0.01 (circles) and At = 0.005 (solid lines). 

the modes n and k are still resonantly coupled, provided 
that no coupling to higher modes exists. However, the 
case of two detuned coupled modes does not apply to the 
scenario discussed here. Decreasing the detuning, i.e. re- 
ducing the value of M, does not only strengthen the cou- 
pling between the modes n = 1 and k = 3 which would 
lead to an exponential growth of the particle number in 
both modes but also enhances the coupling strength to 
higher modes k — 5,7,... because the frequency spec- 
trum becomes equidistant as Af ^ (cf, e.g., [la. l53l|). 
The convergence of the numerical results towards the an- 
alytical expressions for the massless case is demonstrated 
below. 

To study in more detail how the number of produced 
particles depends on the mass we performed numerical 
simulations for a wide range of values for M. The results 
are summarized in Fig. |S1 in a "mass spectrum" where 
the number of particles created in the resonant mode 
A^i (t — 2000) is plotted as a function of M and compared 
to the analytical prediction Eq. H33|) . Particular values 
of M for which Eq. (|32|l gives integer solutions, i.e. exact 
(un-detuned) intermode coupling, are marked by arrows 
and the values of M are indicated. Numerical results 
for these values are not included in the spectrum. Cases 
with exact coupling will be discussed later on. 

The numerical values for iVi perfectly agree with the 
analytical prediction H33|l for values of M larger than 
roughly M = 0.6. For masses smaller than this thresh- 
old value the number of created particles is smaller com- 
pared to the analytical prediction. The mass spectrum 
exhibits a maximum at around M ~ 0.4, i.e. particle 
production in the resonant mode is most efficient for this 
particular mass. When AI < 0.4 the number of created 
particles drops down and approaches the M = result. 
The appearance of a maximum in the mass spectrum 
is clear from the above discussion. For the particular 
value M = 0.4 the equation Sflf = fi^ leads to a value 



150 




FIG. 5; (color online) Number of particles created in the res- 
onance mode n = 1 at time t — 2000 as a function of the mass 
parameter M. The solid line shows the analytical prediction 
Eq. I33II . Arrows pointing towards particular mass values of 
M mark masses for which Eq. 1)33^ is not valid because of 
exact intermode coupling. The coupled modes are given in 
brackets [(1, k)]. No numerical results are shown in the plot 
for those cases. Most of the numerical results are shown for 
different values of the cut-off fcmax to underline stability. 




FIG. 6: (color online) Particle spectra for mass parameter 
M = 0.4 at times t = 1500, 3000, 4500 and 6700. Each spec- 
trum is shown for values fcmax = 10 (dots) and femax ~ 20 
(squares) to indicate numerical stability. 



k = 3.02 which is close enough to the integer solution 
/c = 3 to couple this mode strongly to the resonant mode 
but on the other hand coupling to higher modes is still 
suppressed. Figure shows the particle spectrum ob- 
tained for Af = 0.4 for different times and in Fig. [7| the 
time evolution of the number of particles created in the 
modes k = 1,2,3 and 4 is plotted. For the even modes 
k = 2 and 4 the same oscillating behavior is observed as 
for M = 0.7 and M = 0.2. The couphng of the mode 
k = 3 to the mode n = 1 results in a damping of the 
resonant mode and consequently the number of particles 
produced in the mode n = 1 is smaller than the value 
predicted by Eq. 
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FIG. 7: (color online) (a) Number of particles created in the 
modes k = 1,2,3 and 4 for the mass parameter M = 0.4 
corresponding to the spectra shown in Fig. |H] Part (b) shows 
N3{t) and part (c) N2{t) in each case for the two resolutions 
At = 0.01 (circles) and At = 0.005 (solid lines). 



For increasing masses larger than M — 0.4 the excita- 
tion of higher modes becomes more and more suppressed 
(cf Fig. I2J). Accordingly the numerical results match the 
analytical expression H33|l predicting that the number of 
created particles decreases with increasing mass. De- 
creasing the mass below M = 0.4 enhances the strength 
of the intermode coupling which results in a damping of 
the resonant mode n = 1. Consequently the number of 
particles produced in the mode n = 1 (and also the total 
particle number) is smaller than predicted analytically. 
When studying the limit M ^ the numerical results 
should converge towards the well known results for the 
massless case where all odd modes are coupled [s^ . 
This is demonstrated in Fig. |S1 where the total particle 
number and the number of particles created in the res- 
onant mode n = 1 are depicted for M = 0.2,0.15,0.1 
and 0.05 up to i = 500 and compared with the analytical 
predictions for M = (see also Figs. 4 (a) and (b) 
of llll). While for M = 0.2 the total particle number 
N{t) is still mainly given by Ni{t) a divergency between 
N{t) and N-i{t) starts to become visible for M = 0.15, 
i.e. the influence of the intermode coupling gains impor- 
tance. For M = 0.1 the numerical results are close to the 
analytical M = 0-results and are practically identical to 
them for M = 0.05. 

We now turn to cases with exact coupling between 
two modes. As already mentioned above exact coupling 
of modes takes place if the conditional equation 132|) 
has integer solutions. In particular, exact coupling be- 
tween two modes n and k occurs if Eq. 1)34(1 is satis- 
fied. In [s^ l the authors derived analytical expressions 
(Eqs. (54) and (55) of 33]) for the case that the TE- 
mode il.'^^ (resonant mode) is coupled to the mode 



(5,1,1) 



I.e. 



3O0 

(1,1,1) 



— r2°g ^ is fulfilled. This par- 



ticular case is equivalent to the coupling of the massive 
modes n — 1 and A: = 5 if M = \/2n {Iq = 1). Fig- 
ure El shows the numerically obtained particle spectrum 




500 
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FIG. 8: (color online) Total particle number A'' (circles) and 
number of particles created in the mode fc = 1 TVi (squares) for 
mass parameters M = 0.2, 0.15, 0.1 and 0.05 together with the 
analytical predictions for the massless case Eq. (6.5) (dashed 
line) and Eq. (6.10) (solid line) of [3 to demonstrate the 
convergence of the solutions towards the M = case (see also 
Fig. 4 of 53]). The cut-off parameter fcmax ~ 30 was used in 
the simulations. 
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FIG. 9: (color online) Particle spectra for cj = 2S7? and mass 
parameter M — \/2tt yielding exact coupling between the 
modes n — 1 and k = 5. Dots correspond to fcmax ~ 10 and 
squares to fcmax = 20. 



at four different times. The cut-off parameter fc,nax — 20 
guarantees stability of the numerical results. The nu- 
merical simulations confirm the prediction that practi- 
cally only the modes n = 1 and k — 5 become excited 
and particles are produced exclusively in the two coupled 
modes. Thereby the rate of particle creation is equal for 
the two modes. In Fig. 1101 we show the numerical results 
for Ni{t) and A^5(i) and compare them with the ana- 
lytical expressions Eq. (54) and Eq. (55) of js^] derived 
via multiple scale analysis (MSA). Whereas the numer- 
ical results agree quite well with the analytical predic- 
tion of 33] for long times, one observes a discrepancy 
between the numerical results and the analytical predic- 
tions for "shorter times" up to t '--^ 3000 {ent = 3tt). For 
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FIG. 10: (color online) Number of particles created in the 
modes n = 1 and k — 5 for oj = and M = \Pi-k cor- 
responding to the particle spectra depicted in Fig. |^ The 
numerical results are compared to the analytical predictions 
Eq. (54) [solid hue] and Eq. (55) [dashed line] of 0. The 
numerical results shown correspond to the cut-off parameter 
fcmax = 20 which guarantees stability. 



long times, the analytical predictions nicely reproduce 
the large scale oscillations in the exponentially increas- 
ing particle numbers. For times up to t ~ 500, the nu- 
merically calculated particle numbers grow with a much 
smaller rate than predicted by Eqs. (54) and (55) of [sSj- 
Furthermore, the analytical expressions predict that for 
"short times" iVi and increase with the same rate 
whereas from the numerical simulations we find that the 
production of particles in the mode fc = 5 sets in after the 
production of particles in the n = 1-mode. Apart from 
the differences for short times the numerical results are 
well described by the analytical predictions of [s^- The 
discrepancy between the analytical predictions and the 
numerical results for short times is due to the fact that 
the MSA analysis in [s^ only considers the resonant cou- 
pled modes, but for short enough times all modes should 
be treated on an equal footing [58|. 

As a second example of exact coupling between two 
modes we show in Fi gs. 1111 and El the numerical results 
obtained for M = y/l/Sir for which the mode n = 1 is 
coupled to the mode fc = 4. 

Let us discuss another case with exact coupling of two 
modes which impressively demonstrates that strong cou- 
pling between modes k and I occurs even if Eq. H32|l is 
satisfied only approximately. For M — ^/bn equation 
(|32|l predicts that the mode n = 1 is exactly coupled to 
the mode k = 7 (i.e. fc = 7 is an integer solution of 
317? = nl). The equation 217? ^ fl^ ~ fl^ is not satisfied 
by an integer / but has the solution I ^ 12.04 which is 
close to the integer I = 12. Thus we can expect a cou- 
pling of the mode fc = 7 to the mode 1 — 12. In addition 
one finds that the equation 217° — fi^j — has solution 
m — 16.96, i.e. m ~ 17, and hence / = 12 is coupled to 
m = 17. In the same way the equation 211" = 17^ — fij'y 
which is solved by j ^ 21.93 leads to a coupling between 



FIG. 11: (color online) Particle spectra for lo — 2Qi and mass 
parameter M = •y/7/87r yielding exact coupling between the 
modes n — 1 and fc = 4. Dots correspond to femax = 10 and 
squares to fcmax = 20. 
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FIG. 12: (color online) Number of particles created in the 
modes n = 1 and fc = 4 for uj = 212? and mass parameter 
M — ^/t/Stt corresponding to the particle spectra depicted 
in Fig. ini 



the modes m = 17 and j = 22. Hence from the numer- 
ical simulations we expect to find a particle spectrum 
showing that particle creation takes place in the modes 
fc = 1, 7, 12, 17 and 22. This is demonstrated in Fig. [H 
where the numerically evaluated particle spectrum is de- 
picted for times t = 500, 1000, 1500 and 2000. The cut-off 
parameter fcmax — 50 ensures numerical stability in the 
integration range considered ^. The number of created 
particles Nk{t) is shown in Fig.^jfor the modes fc = 1, 7 
and 17. 

Without having done a detailed analysis we find, as 
a reasonable approximation, that a mode I is (strongly) 



^ From Fig. El 

one observes that also the mode I 
coupled. The equation 2D,^ = D,^ — D,22 t^*^ 
which explains the excitation of the mode I = 27 



27 is weakly 
ion / ~ 26.92 
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FIG. 13: (color online) Particle spectra for ui = 2fii and 
mass parameter M — \/Z-k. Dots correspond to fcmax = 40 
and squares to fcmax = 50. 
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FIG. 14: (color online) Number of particles created in the 
modes n = 1, fc = 7 and I — 17 for lu — 2f2i and mass 
parameter M = ^/En corresponding to the particle spectra 
depicted in Fig. 1131 




M 

FIG. 15: (color online) Number of particles created in the 
resonance mode n = 2 as function of the mass parameter 
M. The solid line corresponds to the analytical prediction 
13311 . Note that the first three values in the spectrum are not 
numerically stable due to an insufficient fcmax. 
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FIG. 16: (color online) Particle spectra for lo — 2Q2 a^nd mass 



parameters M = ^/Ts/Stv and M = 




coupled to a given mode k whenever the ratio |/ — l\/l 
with I denoting the solution of 217° = |r2? ± O^j is of 
the order of or smaller than 10~^, i.e. of the order of or 
smaller than e used in the simulations. 



C. Higher resonance u) = 2fl2 

Now we briefly discuss results obtained for the cavity 
frequency uj — 2fl2- I^i Fig. 1151 we show a numerically 
calculated mass spectrum similar to the one depicted in 
Fig. El The qualitative behavior is the same as discussed 
for the main resonance case. As in Fig. |31 results for 
values of the mass parameter M for which modes are 
exactly coupled are not included in the spectrum but 
marked by arrows with the corresponding coupled modes 
given in brackets. The numerical results again perfectly 
agree with the analytical prediction H33|) for values of 



M larger than a threshold value which is roughly 1.3. 
The maximum in the particle spectrum appears now for 
M ~ 0.8 and the interpretation of the shape of the mass 
spectrum is equivalent to the one given for the case uj — 

In Figs. 1161 and 1171 we finally show numerical results 
for the two mass parameters M = -^/IS/Stt and M — 
y'7/27r yielding exact coupling between the modes 2 and 
7, respectively, 2 and 8. 



VI. PHOTON CREATION IN A 
THREE-DIMENSIONAL CAVITY 

The analogy between massive scalar particles and 
transverse electric photons in a three-dimensional rect- 
angular cavity outlined in Section II allows to interpret 
the presented numerical results as follows: Consider a 
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FIG. 17: (color online) Number of particles created in the 
modes k — 2 and 7 for M — ^/l3/8n, respectively, k — 2 
and 8 for M = \/7/2tt, corresponding to the particle spectra 
shown in Fig. 1161 

three-dimensional rectangular cavity with equally sized 
non-dynamical dimensions ly = Iz = We parametrize 
the size of ^|| in terms of the initial size of the dynamical 
dimension = lx(0) by introducing £ — l\\/lo. If we re- 
strict ourselves for simplicity to the case Uy — Uz = 
the dimensionless mass parameter M reads 

M = /ofc|| = V2(^) . (36) 

Therefore, for fixed any value of M corresponds to 
a particular realization, i.e. size £, of the non-dynamical 
cavity dimensions. We have found that for a particular 
value M the production of massive scalar particles in the 
resonant mode is maximal. Consequently it is possible 
to maximize the production of TE-photons in a three- 
dimensional rectangular cavity by tuning the size t of 
the non-dynamical cavity dimensions. 

For instance, for oj = 2^l\ the creation of massive scalar 
particles in the resonant mode n — 1 la most efficient 
for M ~ 0.4 (cf. Fig. [S]). This corresponds to the 
three-dimensional case with lo — 2^1^-^ ^ and £ ^ 11. 
Hence by designing the three-dimensional cavity such 
that /|| ~ 11^0 the production of TE-mode photons in 
the resonant mode (1, 1, 1) can be maximized. In order to 
maximize the creation of TE-photons in the mode (1, 2, 2) 
when uj — 20^^^ ^ 2) ^^^^ 'll '^^ non-dynamical di- 
mensions has to be doubled, i.e /|| ^ 22 Iq. For uj = 2il2 
we have found that the maximum in the mass spectrum 
is at M 0.8 [cf Fig. 15]. Accordingly, the production of 
TE-photons of frequency fl'^^ 1 1) under resonance condi- 
tions is maximal in a cavity of dimensions /|| ~ 5.6 /o- The 
strong-coupling case ui = 2fli with M — 0.2 where the 
analytical prediction H33|l docs not describe the numeri- 
cal results due to enhanced intermode coupling [cf Figs. 
Hand 121] corresponds to the lowest TE-mode (1, 1, 1) in 
a cavity of size /|| ~ 22 Iq. 

Similarly one can arrange the size of the cavity such 



that particular modes are exactly coupled, i.e. Eq. H34() 
is satisfied. For instance, resonant coupling of the TE- 
modes (1, 1, 1) and (4, 1, 1) corresponds to w = 217° with 
M — -y/T/STT [cf Figs. 11 and 12] and is therefore realized 
in a cavity of size /y ~ 1.5 Zq. Finally, choosing Z|| ~ 

0.63 i-e. M = VSvr, couples the TE-modes (1,1,1), 
(7, 1, 1),(12, 1, 1), (17, 1, 1) and (22, 1, 1) in the resonance 
case UJ = 2^l^^ ^ [cf Figs. 13 and 14]. 

In summary, the mass spectrum Fig. \^ can be in- 
terpreted in the following way if we set n|| = 1 and 
Lu — 2n^^ ^ j^y. For Z|| = Zq, i.e. cubic cavity, the modes 

(1,1,1) and (5,1,1) are resonantly coupled (cf. Figs. 
and I10|l . Enlarging /|| with respect to Iq increases 
the production of resonance mode photons (1, 1, 1) until 
Z|| = 1.5 Iq (M — y/7/8n) is approached where the modes 
(1, 1, 1) and (4, 1, 1) are exactly coupled (cf. Figs. [TTland 
I12|l . When increasing further photon creation in the 
TE-mode (1, 1, 1) becomes more and more efficient and 
is perfectly described by Eq. Ip!^ . Reaching Zy ~ 7.4 Zq 
(M ~ 0.6, the threshold) the intermode coupling starts to 
become noticeable causing slight deviations of the numer- 
ical results from the analytical prediction. For Zy 11 Iq 
(M ~ 0.4) the production of TE-mode photons is most 
efficient. The number of Photons created in the mode 
(1,1,1) is smaller than the analytical prediction Eq. I|33() 
because of the couphng of the modes (1,1,1) and (3,1,1). 
When increasing ?|| beyond ~ IIZq the strength of the 
intermode coupling is enhanced drastically and conse- 
quently the number of produced TE-mode photons de- 
creases rapidly. For Zy ~ 22^0, for instance, the mode 
(1,1,1) is (strongly) coupled to the modes (3,1,1) and 
(5, 1, 1) (cf Fig. 121). Reducing Z|| with respect to Iq (i.e. 
going to masses M > V27r) lowers the efficiency of pho- 
ton creation in the resonant mode. The mass spectrum 
Fig. 1151 owns an equivalent interpretation. 



VII. CONCLUSIONS 

The production of massive scalar particles in a one- 
dimensional cavity, or analogously the creation of TE- 
mode photons in a three-dimensional rectangular cavity, 
has been studied numerically for resonant wall oscilla- 
tions. 

We have found perfect agreement between the numer- 
ical results and analytical predictions of js^ in the case 
that no modes are (strongly) coupled. When two modes 
are exactly coupled, i.e. Eq. (|32|l possesses integer solu- 
tions for I and k, the numerical results agree with ana- 
lytical predictions of for sufficiently long times but 
disagree for short times. The discrepancy for short times 
is ascribed to properties of the multiple scale analysis 
used in [ssl] . 

The effect of the intermode coupling has been studied 
in detail which is only possible by means of numerical 
simulations. As main result we have found that a particu- 
lar mass exists for which the production of massive scalar 
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2fi°^ ^ has been dis- 



particles is most efficient. The appearance of a maximum 
in the mass spectrum, i.e. the number of created parti- 
cles after a given time as function of mass, is explained by 
the increasing strength of the intermode coupling when 
decreasing the mass below a certain threshold value. 

The analogy between massive scalar particles in a one- 
dimensional cavity and TE-mode photons in a three- 
dimensional cavity allows the conclusion that the effi- 
ciency of TE-mode photon production from vacuum in a 
resonantly vibrating rectangular cavity can be controlled 
(and maximized) by tuning the size of the cavity when 
keeping the quantum numbers ny,nz corresponding to 
the non-dynamical cavity dimensions fixed 

The main resonance case lu 
cussed for a cavity with equally sized non-dynamical di 
mensions — ly — Iz in detail in section VI. We have 
shown that photon creation in the resonant mode (1,1,1) 
is most efficient if the size l\\ of the non-dynamical cav- 
ity dimensions is ~ 11 times larger than the dynamical 
cavity dimension. The existence of a certain cavity size 
which maximizes photon creation is among other things 
explained by the fact that intermode coupling takes place 
even if Eq. H32I) is satisfied only approximately. If l\\ is 
larger than this value the intermode coupling is so strong 
that higher frequency modes hke (3, 1, 1) and (5, 1, 1) cou- 
ple to the resonant mode (1, 1, 1) and strongly damp its 
evolution. Furthermore the coupling of particular field 
modes by tuning the size of the cavity has been studied. 
The effects provoked by the intermode coupling can be 
studied in full detail only by means of numerical meth- 
ods. 

Our findings demonstrate that the intermode coupling 
in dynamical cavities plays an important role. Even if 
analytical results are known, numerical simulations are a 
very useful and indeed necessary tool because only they 
can completely take into account the intermode coupling. 
In order to study photon creation associated with the full 
electromagnetic field in a dynamical cavity also the con- 
tribution of the transverse magnetic modes (TM) to the 
photon production has to be considered. Studying TM- 
modes numerically represents a more demanding task be- 
cause of the more complicated so-called generalized Neu- 
mann boundary condition these modes are subject to. 
This will be addressed in a future work. 
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APPENDIX A: REMARKS ON NUMERICS 

To solve the system of differential equations formed by 
Eqs. H22|l and H23|) numerically we decompose S,^\t) and 
"nrT^t) in their real and imaginary parts: 



e 



(m) 



(m) (m) 
n t In 



(Al) 



The resulting coupled system of first-order differential 
equations can then be written in the form 



(A2) 



with real vectors X}'^\t) and matrix W_{t). Choosing 
the representation 

vim) _ ( {m) (m) (m) (m) (m) (m) (m) {ra)-.T 

.(A3) 

where we have truncated the infinite system via intro- 
ducing the cut-off parameter K = fcmax, the x AK - 
matrix W_{t) becomes 



W(i) = - 



C~{t) C+(t) ~A+{t) A-{t) 

C+{t) C-{t) ~A-{t) A+{t) 

A+{t) -A-{t) C-{t) C+{t) 

lA-{t) -A+{t) C+{t) C-{t) 



(A4) 



with the K X K - matrices C^(i) = [cf„{t)] , 1 < fc,n < 
K and diagonal matrices A'^{t) = [a„„(i)] where a^„(t) 
and c^i^{t) are defined in Eq. TI^ and Eq. (|2SJl, respec- 
tively. The number of particles (|18|l created in a mode 
n at t = ti may now be expressed in terms of the real 
functions: 



m=l ™ 

+ [A-{t^)v(rHh) + A+{h)i:"Hh)f} (as) 



which in the particular case ti — AfT with T the period 
of the cavity oscillations and integer Af reduces to 

m—l 

(A6) 

In order to calculate (|A5(I the system (|A2ll has to be 
evolved numerically if-times (m is running from 1 to 
K = /cmax) up to t = ti with initial conditions 



z;(")(0)=x(™)(0) = y(™)(0) = 



(A7) 
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FIG. 18: (color online) The function dk{t) [Eq. lIXTTl l for 
LO = 20? and M = 727r [panel (a) and (b)] and M = 0.4 
[panel (c) and (d)]. In any case dk{t) is shown for k = 1, 5 
(upper bands) and the last five values k — fe^ax — 4, fcmax 
(lower bands). With "err" we denote the presetted values for 
the relative and absolute error used in the numerical simula- 
tions performed with, in these cases, the Runge-Kutta Prince- 
Dormand method. 



and 



(A8) 



Besides investigating the stability of the numerical solu- 
tions in dependence on the cut-off K the quality of the 
numerical solutions can be assessed by checking the va- 
lidity of the Bogoliubov relations 

V [A,Uti)AUiti) - B*,^{h)B„,k{ti)] = Snk (A9) 



m 



[A„Uti)B*^kiti) - B*^^{h)Amk{ti)] = 0. (AlO) 



In order to solve the system ()A2|I numerically we applied 
integration routines based on different standard solvers. 
Mainly employed were the Runge-Kutta-Fehlberg 4th- 
5th order method (rkf 45) and the Runge-Kutta Prince- 
Dormand method (rkSpd). Source codes provided by the 



GNU Scientific Library (GSL) Ha as well as the MAT- 
PACK - Library 6Qj were used. 

In the following we discuss the accuracy of the numer- 
ical simulations by considering the quantity 



dk{t) = l-J2{\-^rnk{t)\^-\B,nk{t)\' 



(All) 



indicating to what accuracy the diagonal part of the rela- 
tion l|A9p is satisfied by the numerical solutions. Generic 
examples for dk{t) are shown in Fig. El Panels (a) and 
(b) correspond to the exact coupling case M — \/2tt (cf. 
Figs. l^and [Tn|l with the absolute and relative errors (err) 
for the Runge-Kutta Prince-Dormand method (rkSpd) 
[59I I have been set to 10~^ (a) and lO^^^ (b) Thereby 
two "bands" are shown. The upper one correspond to 
fc = 1 to 5 whereas the lower one correspond to fc = 16 to 
^max = 20. The deviation from zero is larger for higher 
k because these modes are more affected by the trun- 
cation of the infinite system through the cut-off fcmax- 
Comparing the absolute value of the maximal deviation 
of dk{t = 8000) from zero which is ~ 3.5 x 10~^ for 
err=10^^ and ^ 1.5 x 10^^ for err=10~^^ with the num- 
ber of particles created in the resonantly excited modes 
Ni{t = 8000) - 350 and N^it = 8000) ~ 400 demon- 
strates that the numerical simulations guarantee a good 
accuracy. 

In panel (c) and (d) of Fig. E| we show dk (t) for the 
case M = 0.4 (cf. Figs. [HI and [7|) for the cut-off values 
fcmax — 30 (c) and fcmax = 50 (d). The numerical sim- 
ulations have been performed with err=10^^ and again 
two bands are shown corresponding to the first five (up- 
per band) and last five (lower band) values of k. For 
fcmax = 50 the deviation of the absolute value of dk from 
zero for the last values k = 46, ...,50 [panel (d)] is slightly 
larger compared to the deviation for the last five modes 
for /cniax = 30. But the deviation of dk{t) from zero for 
the first modes k — 1, ...,5 is smaller for /cmax = 50 than 
for /cmax — 30, i.e. the accuracy for the first modes im- 
proves when increasing femax as it is expected. Comparing 
\dk{t = 2000)1 - 3 X lO"'' for fc^ax = 50 with the number 
of created particles Ni{t = 2000) ~ 124 demonstrates 
again the accuracy of the numerical simulations. The re- 
maining Bogoliubov relations are satisfied with the same 
accuracy. 
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